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Abstract: An exact solution is presented for a swinging Atwood's machine. 
This teardrop-heart orbit is constructed using Hamilton- Jacobi theory. The 
example nicely illustrates the utility of the Hamilton- Jacobi method for finding 
solutions to nonlinear mechanical systems when more elementary techniques 
fail. 
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A swinging Atwood's machine (SAM) is an ordinary Atwood's machine in 
which, however, one of the weights can swing in a plane [1]. In the polar 
coordinates (r, 9) , the nonlinear Hamiltonian describing the dynamics of SAM 
is 

H = T + V = -(m + M)f 2 + -mr 2 9 2 + gr(M + mcos9), (1) 

where M is the mass of the nonswinging bob, m is the mass of the swinging 
bob, g is the gravitational acceleration, and the polar angle 9 is measured from 
the negative vertical y-axis as shown in Figure 1(a). When the mass ratio 
/i = M/m = 3, numerical studies of SAM reveal an unexpected feature. To wit, 
all orbits which start at the origin (with initial conditions r = e, f = v, 9 = 
9o, — 0) execute one symmetrical loop and return to the origin, no matter 
what the launch angle (9 ) or speed (v). When the launch angle is small, the 
orbit is called a teardrop, and when the launch angle is larger (say 9$ > tt/2) 
the orbit is called a heart (see Figure 1) [2]. 

Since the SAM Hamiltonian is nonlinear, we generally expect that it will be 
impossible to find (simple) exact solutions for SAM [3]. In this paper I show 
that the "teardrop-heart" orbit of SAM is an exception to this rule. Specifically, 
I use Hamilton- Jacobi theory to construct an exact solution for the teardrop- 
heart orbit. This example would make a nice addition to an advanced mechanics 
class since it illustrates the power of the Hamilton- Jacobi method to arrive at 
an exact solution which apparently alludes more elementary methods. 

Recall that the game of Hamilton-Jacobi theory is to find some canonical 
transformation which will separate variables in the Hamilton-Jacobi equation. 
A fortuitous choice for the SAM Hamiltonian is the point transformation 

r = \(e+V 2 ), (2) 

9 = 2arctan[(£ 2 - r? 2 )/2^]. (3) 

This transformation compresses the angular variable by a factor of two, and 
then changes from polar coordinates to parabolic coordinates with parabolas 
centered about the horizontal x-axis [4] . The inverse transformation is 



£ = Wl + sin(0/2), (4) 

r? = ±Wl -sin(0/2). (5) 

If I further set M — 3m (and M — 1) as is done throughout the rest of this 
paper I find that the generalized momentum in the new coordinate system is 

P 5 =m 2 + V 2 ) = ^, Pn=m 2 + V 2 ) = ^, (6) 
from which the (time-dependent) Hamiltion- Jacobi equation is calculated as 

-2 5 (£ 4 + 77 4 )} =0. (7) 



as i I i 

dt + (e+V 2 ) I 8 



dsy /as 
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The equation for the generating function is now separable and can be ex- 
pressed in the form 



(8) 



Separating the time and space parts of Eq. (7) I find that the first separation 
constant is just the total energy, E: 



dS 
~dt 



-E, 



OS 



(e+v 2 ) \8 

Integration of Eq. (9) yields 



fdSY 



S(t) = -Et. 



Equation (10) can be written as 



as 
at 



asy 



(9) 
(10) 

(11) 
(12) 



Next I can separate the spatial part of the Hamilton-Jacobi (Eq. (12)) equation 
giving 



as 



m g e - m 2 = i, 



(13) 



(14) 



where I is the spatial separation constant. The complete solution to the Hamilton- 
Jacobi equation is then given by 

S = J(I + 8E£ 2 - 16 3 £ 4 )^ + [(-! + 8EV 2 - 16gv*)*dri - Et. (15) 

To "solve" the original problem I do not need to integrate Eq. (15) directly, 
but rather the orbit equation: 



as _ i r i 

97 ~ 2 J (I + SB!? - IQg^i) 



s 2 J (_/ + 8Er] 2 - Wgrj 



-dr/ = a (16) 



which provides an implicit relation between £ and rj (and, therefore, r and 9), 
and the parametrized trajectory equation: 



-s 



4£ 2 



as 



dE J (I + 8EC 2 - 16 5 £ 4 )5 



4r/ 2 



(-I + SErj 2 - 16 OT 4 )2 



T dri = t + (17) 
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where a and (3 are constants of integration which may, in general, be a compli- 
cated mix of the original initial conditions. 

Just as important, I can also solve for the new constant of the motion, /, by 
adding Eqs. (13) and (14), and then eliminating the energy constant to get 



/ = 2(e + v 2 )(v 2 e eV) + wwue v 2 )/(e + v% 



(18) 



In the original polar coordinate system (apply Eqs. (4) and (5) to Eq. (18)) the 
new constant of motion reads 



I(r,f,6,6) = r 2 9 



r9 

r,os,-l- T sinl- 



gr 2 sin ( \ cos 2 I - 



(19) 



In principal the orbit and trajectory equations arising from the generating 
function represent a complete solution to the original problem — in this case 
finding the motion of SAM when \i = 3. In practice, though, there are still two 
difficult obstacles to overcome. First, actually finding explicit integrals to the 
orbit and trajectory equations may not be possible, and second, the resulting 
relations are often, at best, only implicit functions which can not be inverted 
for the relevant variables. 

In this particular example, the first of these obstacles can be overcome — 
the orbit and trajectory equations can be explicitly integrated with Elliptic 
functions. However, instead of using special functions to arrive at a solution, I 
would next like to turn my attention to a specific solution — the teardrop-heart 
orbits — and show that for this special class of orbits it is possible to obtain a 
simple and elementary solution for the motion of SAM. 

The key observation used in arriving at this elementary solution is as follows: 
for all teardrop-heart orbits the constant of motion "I " is equal to zero. To see 
why this last observation is true consider the teardrop-heart orbit as it crosses 
the negative vertical axis (see Figure 1). At the axis 9 = and r = 0. Plugging 
these values into Eq. (19) it follows that I(r,r,9,9) = for all t. Next I show 
that when 1 = the integrals in the orbit and trajectory equations can be solved 
by elementary functions, and the resulting functions can be inverted to arrive 
at a simple orbit equation relating r and 9. 

Examining the special case / = 0, the orbit equation (Eq. (16)) becomes 



2a 



= f 



1 



(8££ 2 - 16 5 £ 4 )2 2 J (se v 2 
which can be integrated straight away to give [5] : 



16 OT 4 ) 



-di] 



2a = 



'8E 



2E - 4g£ 2 



'2E 



HO 



In - Agrf 



(20) 



*2e) - lnfa)] } 
(21) 
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After some additional algebraic manipulations I find the implicit orbit equa- 
tion in £ and w to be 



where the positive constants a and k are defined by 



(22) 



a = e 2a ^, k = % (23) 

Using the inverse point transformation Eqs. (4) and (5) applied to Eq. (22) I 
get the implicit orbit equation in polar coordinates as 



v 1 ~ 8in 2 v/l-fcr(l + shig)-l ^ 



^1 + sinf y/\-kv{\-save)-\ 



At this point it is helpful to pause and check this result by considering the 
limiting case of a teardrop orbit. Namely, the solution of Eq. (24) in the limit 
of small 9. Expanding Eq. (24) in 8 about zero I find [6]: 



±a = 1 + ^T=Tr +0{e2) - <25) 



Or, written another way 



r(9) « r 



2 

/ li \ 

1 - 



r Q = 1/fc, d Q = 2(-l±o), (26) 



which indeed produces the teardrop shaped orbit shown in Figure 1(a). 

Now I return to the main stream of my presentation by attempting to 
find an explicit orbit equation. Rewrite Eq. (22) in the form a£yl — krj 2 — 
r]\/l — fc£ 2 = a£ — T], square both sides, and collect like terms to discover 
2a — k£r/ — ka 2 ^ri = 2a\f \ — fc£ 2 — krf + k 2 ^ 2 r\ 2 . Next notice (from Eqs. (4) 
and (5)) that £77 = r cosfl. Evidently, (fc 2 + k 2 a A - 2a 2 k 2 )r 2 cos 2 (§) + {-Aak - 
4a 3 fc)r cos (|) + (4a 2 fc)r [l + sin (§)] + (4a 2 fc)r [l - sin (|)] = 0. Additional al- 
gebraic simplification results in the exact explicit teardrop -heart orbit equation 

4a (a 2 + 1) Icos (|)| - 2a 
k (a 2 - l) 2 cos 2 (I) 

that is, the teardrop-heart orbit occurs for positive values of r, and this condition 
places the above restriction on 9 where (at r = 0) 

9 a = 2 arccos ( ^ ^ - j . (28) 
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Also, when the orbit crosses the negative vertical axis its length is (at 8 = 0) 

4a 1 

r ° = T ' (^TIF (29) 

That's it. 
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and 




1 + 







+ 



(-2 + 2kr-2Vl-kr + krVl-kr) 2 



+ 0(9 3 ). 



2y/l~kr 



SVl-kr (-1 + fcr) (l + Vl-fcr) 
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Figure Captions 

Figure 1: (a) Teardrop, (b), (c) heart and (d) a rabbit ear orbit of a swinging 
Atwood's machine (SAM) when \i = M/m = 3. 
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